Multipartite continuous- variable entanglement distillation using local squeezing and 

only one photon-subtraction operation 
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In this paper, we study entanglement distillation of multipartite continuous-variable Gaussian 
entangled states. Following Opatrny et al.'s photon subtraction (PS) scheme, the probability of 
successful distillation decreases exponentially with the number of parties N. However, here, we shall 
propose an entanglement distillation scheme whose success probability scales as a constant with TV. 
Our protocol employs several local squeezers, but it requires only a single PS operation. Using 
the logarithmic negativity as a measure of entanglement, we find that both the success probability 
and the distilled entanglement can be improved at the same time. Moreover, an TV-mode transfer 
theorem (transferring states from phase space to Hilbert space) is presented. 

PACS numbers: 03.67.Mn, 03.67.Hk, 42.50.Dv 
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Entanglement, particularly multipartite entanglement 
state, is one of the most fundamental and puzzling aspect 
in quantum mechanics. However, entanglement is such 
a fragile resource that it may be easily degraded dur- 
ing its interaction with the environmental noise. To this 
point, entanglement distillation (always in a probabilistic 
way) has been proposed to increase the entanglement in 
the noise-disturbed entangled state [1, 2]. Restricted by 
the famous No-Go theorem in continuous variable (CV) 
entanglement distillation[3~5], lots of efforts have been 
devoted to the non-Gaussian operations. As an exam- 
ple, photon subtraction(PS) operation[6], proposed by 
Opatrny et al in 2000, is principally simple and can be 
readily implemented with beamsplitter and photon detec- 
tors. Very recently, about 10 years after Opatrny et aVs 
pioneering work, an experiment which faithfully imple- 
ments the PS-based two mode entanglement distillation 
has been reported [7] . One of the challenge in this experi- 
ment is the extremely low successful probability, which is 
mainly due to the rathcr-high-transmittance beamsplit- 
ter used in PS operation — For one thing, the beamsplit- 
ter must own a relatively high transmittancc to guarantee 
an entanglement-enhanced distillation [8]. For the other, 
high transmittance means low reflectiveness in beamsplit- 
ter and hence, low probability in PS operation. Assum- 
ing the beamsplitter's transmittance is 0.90, as shown in 
Rcf.[7], detectors 's efficiency 0.10, the probability of each 
successful local PS operation is upper-bounded by 10~ 2 . 
The successful probability of the whole entanglement dis- 
tillation will be even lower, decreasing exponentially with 
the number TV of local PS operations: 10~ 2N . This is re- 
ally a serious problem if TV-partite continuous variable 
entanglement distillation is involved. 

In this paper, we consider the distillation of TV-partite 
CV gaussian entanglement state with only one-time pho- 
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ton subtraction. Generally for TV partite (especially, the 
symmetric) Gaussian state, one-time PS will modify the 
permutation symmetry and the distilled entanglement, 
measured with logarithmic negativity (log-neg), will be 
even worse. However, if assisted with local squeezing, we 
show that the entanglement can be improved [9]. More- 
over, one-time PS can give a substantial increment in 
distilling probability, which keeps constants O(10 -2 ) for 
arbitrary partite number TV. 

Our entanglement distilling scheme is briefly shown in 
Fig.l, where Fig. 1(a) describes the typical TV-time PS 
scheme which is a straightforward generalization of Opa- 
trny et al's scheme to TV partite case. Fig. 1(b) is the our 
new distilling protocol with local squeezing ( described 
by symplectic operation S(si) ) and only one-time PS. 
Throughout this paper, for convenience, we assume that 
the second input-mode Bi of beam-splitter are vacuum 
modes, representing a simple and concise expression of 
PS operation in Phase space(See AppendixA). 

Preliminaries. Our result can be conveniently derived 
in Phase space. Let's now introduce the basic facts 
and notation of CV TV modes state. First, it is con- 
venient to express each mode, say k mode, with the 
field quadrature operators £k = (a* + a£)/\/2,pfc = 

(afc — Oj.) / (iy/2) , with a/b, a k being the mode annihilation 
and creation operators. By defining the vector of quadra- 
ture operators X = (xi,pi,--- ,xn,pn), the commuta- 
tion relation can be written as [X m ,X n ] 

" = J' 
TV-mode system p which resides in infinite-dimensional 
Hilbert space can be conveniently represented by the 
characteristic function in 2TV-dimensional real vector 
space, i.e., phase space x(0 = Tr[pexp(iA" T £)], £ £ 
R 2N . Gaussian state is a special kind of quantum state 
whose characteristic function is Gaussian in phase space, 

inwhich x is 



ifi m „, with 



The density matrix of arbitrary 



X (0 = exp[-±£ T V£ + zx T £] ,£ e 
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FIG. 1. Multi-partite Entanglement distillation scheme with Photon Subtraction technique. (a)Typical PS-based TV-partite 
distillation with TV times photon subtraction. The successful distillation is heralded by the events when all the detectors 
register non-vacuum result, with the probability scaling as 10~ 2Ar . (b) One-time PS and local squeezing assisted entanglement 
distillation. The success probability conditioned on only-one detector's detection result keep fixed in O(10~ 2 ). 



quadrature average x = Tr[TVp] and V denotes the co- 
variance matrix V; m = \(XiX m + X m X t ) - {X l )(X m ). 
In the following, we also use the wigner function, the 
Fourier transform of x(0 to express the state evolu- 
tion. A standard normalized wigner function is defined 
by W(r) = J R2N (g^exp [-ir T S;] X (£), with n e K 2 

andr= (n,f 2 ,--- ,r N ) eR 2N . 

Multi-partite entanglement states. Multipartite CV 
Gaussian entanglement plays a prominent role in fu- 
ture quantum network and quantum communication 
protocol[10, 11]. In this paper, we are mainly interested 
in the distilling of a family of genuinely N partite sym- 
metric Gaussian entangled state[12-14]. The covariancc 
matrix is given by 



V N =e® \T){1\ + (a-e)<g> J 



AT, 



(1) 



inwhich \T) is the unnormalized TV dimensional real vec- 
tor \T) = (1,1,1- •• ,1) T e R", and I N is the identity 
matrix in TV dimensions, a = diag(a 1 b),e = diag(c,d) 
are 2x2 diagonal matrix 

0= — (e 2ri + (TV - l)e~ 2r A , c = — (e 2 '' 1 - e^), 
2TV v v ; ' 2N ' 
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The genuinely TV partite state in Eq.(l) can be ex- 
perimentally prepared with a particular sequence of 
TV — 1 phase-free beam splitters and TV squeezed input 
states[12]. To be simple, we will mainly focused on 
the unbiased states, namely n = ri(r 2 ) = \ ln[(TV — 
l)sinh(2r 2 )] + ±ln[Vl + P - 1) sinh(2r 2 )]- 2 + 1]. 



Distillation with local squeezing and one-time PS. Let's 
now derive the state evolution of the our one-time 
PS distilling protocol. As shown in Fig. 1(b), N local 
squeezing symplcctic transformations S{si)i = i l ... n arc- 
applied before PS. This corresponds in phase space to 
a transformation of covariance matrix Vn ~ > = 
[®? =1 S{s. l )}V N [®fL 1 S(s i ) T }. Without loss of general- 
ity, we assume that only PS operation is performed 
in the A\ mode. Implemented with a beam-splitter 
(transmittance T), the PS operation (Fig. 1(b) inset), 
couples the TV mode Gaussain state V$ with the vac- 



uum mode B\. 

V BS = B(Vg ( 



matrix B 



Tin 



The (TV + 1) mode states now follows 
\h)B T , with B being the symplcctic 

-Vri 2 \ 

^2(jv-i) • Finally, a suc- 

'RI 2 VTl 2 J 

ccssful distillation is heralded if the detector register non- 
vacuum results. According to Appendix A, one can find 
the distilled state is a linear combination of two gaussian 
state 

Pdis = T—rP( T i) ~ T—TP( r 2): ( 3 ) 



V 



V 



with 5 = y / det(r 2 + 7 2 /2) and p(T) being a normal- 
ized TV-paritie gaussian state with covariance matrix Y . 
The Ti and T 2 arc defined by partitioning of matrix 

Vbs = ( ^ T A l ) ,r 2 = ri - M(A + 7 2 /2)- 1 M T , 

where T X ,M, A are 2 TV x 2 TV, 2TV x 2, 2 x 2 matrice re- 
spectively. The success probability of distilling follows 
P succ = (5 — 1)/ 5. With the transfer theorem form phase 
space to Hilbcrt space (See Appendix B), one can easily 
compare the entanglement before and after distillation. 



3 




0.05 0.1 0.05 0.1 

r r 



FIG. 2. Entanglement before and after distillation for N = 3 
partite genuine Gaussian entanglement state (Eq.(l)) with 
T2 = r — 0.05, n = ri(r 2 ) = 0.099 and si = S2 = S3 = 
s,T — 0.9.(a)Entanglement distillation of unbiased entangled 
state. There exists an optimal local squeezing s which opti- 
mizes the distilled entanglement. For r — 0.05, s opt = 0.07 
and corresponding maximal log-neg and success probability 
is 0.9818 and 9.67 x 10~ 3 . (b) The probability of successful 
distillation as a function of squeezing s. (c)Optimized local 
squeezing s opt as a function of initial squeezing parameter r. 
Sopt scales linearly with the initial parameter r: s op t ~ lAr. 
(d) The corresponding entanglement with s op t- In the nu- 
merical simulation, we truncate the photon number of each 
mode at D = 7, namely, we consider only the contribution of 
|0), |1), 1 2), ■ ■ ■ |6). By choosing D = 7, the error during our 
state transfer can be controlled with in 5 x 10 -5 . 



In Fig. 2, we evaluate the entanglement after and be- 
fore distillation with the log-neg [15] as the figure of merit 
for entanglement. For simplicity, we consider three-mode 
entanglement distillation as the example. For a fixed ini- 
tial squeezing r = 0.05, as shown in Fig. 2(a), we plot 
the entanglement as a function of s\ = S2 = S3 = s. 
Through our simulation, we assume the transmittancc 
of PS beamsplitter is T = 0.90 which is available in re- 
cent experiments [7]. The probability of success is briefly 
shown in Fig. 2(b). It should be noted that the success 
probability which is about O(10~ 3 ). This is mainly due 
to the rather low initial squeezing (r = 0.05), which re- 
sults extremely low photon number in each transmission 
mode, which certainly decrease the probability of being 
photon subtracted. Our method can be applied for even 
stronger squeezing. In Fig.(c), we increase r and find 
the optimal squeezing s op t which may maximize the log- 
neg of output entanglement state. The numerical results 
support the linearity reliance of s op t upon the increas- 
ing r. Also, we plot the corresponding optimized log- 
neg and success probability in Fig. 2(d) and Fig. 3. The 
success probability (Fig. 3) is about 0(1O -3 ) which is an 
pronounced improvement compared with the 3-time PS 
stratcgy(O(10- 9 )). 

Discussions and Prospectives. We presented here a 




FIG. 3. All parameters are choosen as in Fig. 2(d). Only here 
we use the optimal local squeezing s opt and plot the success 
probability. 



photon-subtraction based entanglement distillation for 
arbitrary N— partite continuous variable entanglement 
state. As an example, now in this paper, only the 
three-partite symmetric Gaussian state is involved. This 
method is applicable for arbitrary N~ partite CV state. 
Indeed, even for N = 2, this improvement in both log- 
neg and success probability also applies. As an auxiliary 
result, we also derive the transfer theorem for N— par- 
tite Gaussian state from Phase space to Hilbert space. 
We can envisage that this theorem could find more ap- 
plication in the entanglement evaluation tasks, such as 
entanglement swapping and entanglement distribution. 
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Appendix A: Wigner function description of Photon 
Subtraction 

In ideal perfect photon subtraction is described 

with the annihilation operation: a = X)»^Li \fn\n — l)(n| 
in Hilbert space. However, this is not an unitary op- 
eration and cannot be implemented deterministically. 
A convenient way is to use beamsplitter and photon 
detectors[6]. For ease, we consider a single partite state 
Pi as input, the PS operation II° n (see Fig.l) can be rep- 
resented with a completely-positive map from pi to nor- 
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malized output state pi = £(pi)/Tr[£(pi)], with 
£{ Pi )= Tr B [U AB { Pi ® |0} Bi (0|) U[ B (l Ai ® fir 



where 11° 



(n\ denotes the positive oper- 



ators projecting to non-vacuum subspace and U AB = 
cxp[arccos(v / T)(a| 1 . a B — a A .a g )] denotes the Beam- 
splitting operation between Ai and Bi modes. 

In our calculation, indeed, it is convenient to use 
the wigner function to describe the PS process above 
(Eq.(Al)). In fact, the operator fl° n = I - |0) B .(0| is 
a difference of two operation whose wigner function are 
both Gaussian[16], i.e., 



W(f h n on ) = — (1 - 2 exp[nl 2 ff ]) (A2) 

In case that the input state pi is Gaussian, the wigner- 
function of distilled entanglement can be easily formu- 
lated with the linear combination of a series of Gaussian 
function, each of which can be conveniently expressed 
with the covariance marices. 



Appendix B: Quantum State from phase-space to 
Hilbert Space 

In this section, we give the detailed techniques we 
use in the processing of multi-partite Gaussian quantum 
state from phase-space to Hilbert Space. For a iV-partite 
CV state, the density matrix follows 



P 



^Tr[pD(fl)]D(-iT),fl= (pi,p 2 ,--- ,P-n) e C N 

7T 



with D{p) = exp [pa^ — p*a\ = exp [(/i, p,*){cft , — a) T ] 
being the N mode displacement operator. 

The matrix entries of p can be conveniently obtained 
by observing the equation[9]: 

(fci,fc 2 ,--- 1 k N \D(-p)\m 1 ,m 2 , ■ ■ ■ ,tojv) 

nti^nti^ 1 



: CXp 



tt' - tp + t'p* 



where t = (ii, t 2 , ■ ■ ■ ,tjv)j f = {^i^ 2 ^ • • • > *jv) is the N- 
dimcnsional real vector. 
By noticing the fact that 



(a\ -a) = iL N (xi,p 1 ,x 2 ,P2, ■ ■ ■ ,XAr,p w ) T 

l n =((^- ^)^i N )p,p kl = { 6 ; k - 1 ' 1 

one obtains that (after integrating the p) 

(kx,k 2 ,--- ,kisr\p\mi,m 2 , ■■■ ,m N ) 



(Bl) 

k < N 
k > N 



n<=iMiisr%i>= 



\t=t'=0 



inwhich 



F= exp 



-(t,t')R(t,t') T 



/Vdet(r + 7 2JV /2) 



(B2) 



(B3) 



R= a x ®I N + (a z ® /jy) i^(r + I 2N /2)- L V N (a z ® /jy) 

(B4) 

and <J Xl a z is the Pauli matrices. Then, one can check 
that the state p(Eq.(B2)) is now automatically normal- 
ized, i.e., Tr[p] = 1. 
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